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Abstract 

In a recent publication it has been shown how to generate derivatives with respect to 
atom coordinates of Slater-Koster matrix elements for the tight binding (TB) modelling of a 
system. For the special case of a mixed second partial derivative at coordinate singularities 
only the results were stated in that publication. In this work, the derivation of these results 
is given in detail. Though it may seem rather 'technical' and only applicable to a very special 
case, atomic configurations where the connecting vector between the two atoms involved in a 
two-centre matrix element is aligned along the z-axis (in the usual approach) require results 
for precisely this case. The expressions derived in this work have been implemented in the 
DINAMO code. 



1 Slater-Koster matrix elements 

In [1] it was shown how the Slater-Koster matrix element between orbitals locahsed at two 
different atoms could be systematically expressed in terms of the distance R between the atoms, 
and the Euler angles a and /? describing the orientation of the connecting vector R = (X, Y, Z) = 
R{cos a sin /3, sin a sin /3, cos f3) in space and parametrising the rotations of the original coordinate 
system into a new one necessary for the evaluation of the matrix element. The general expression 
for a matrix element {hrriil H\l2m2), according to [1], with /< = min(/i,/2), is 

l< 

{limi\H\l2m2){a,(3,R) = ^ 

m'=l 

X {hl2\m'\){R) + 2^^,(a)yl^,(a)d|;^^|o(/3)d|^^|o(/3)(/i/20)(i?) 

Here Ihrni) denotes a state whose wavefunction has an angular dependence characterised by a 
real spherical harmonic Yi-m^ on atoms i = 1,2, respectively. The real spherical harmonics are 
defined as 

Yim = S^oYio + (1 - 5^o)V2i-ir [r{m)ReYnm\ + r(-m)Imyi|^|] (2) 

with r(m) = 1 if m > and r(m) = if m < 0. The Yim are the ordinary complex valued 
spherical harmonics, with the phase convention Y^'^{9,Lp) = (— l)'"y;_m(0, (/?). 
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in is an operator that is cylindrically symmetric about all possible orientations of the 
connecting vector R; it can be the Hamiltonian, but also the identity, in which case the matrix 
elements would be overlap integrals between localised states. The other quantities appearing 
are 

I (— 1)™ rr(m) cos(|m|a) + r(— m) sin(|m|a)l if m 7^ 
^-(a):=|_i^ if^ = 

„ . . f (—I)™' Trf—m) cosdmla) — rfm) sin(|m|a)l if m 7^ 
BM--=[l ^ L I U .^^^^ (4) 

^rnm' '= ^m[( — 1)™" 4m|m' + ^|m|-m'] (^) 



\m\m' ' \m\—in'l 

/' - ^' 1 

''|m|m' "Iml-m'J 



and the Wigner d-function [3] 



= 2-'(-l)'-™' [(/ + m)!(/ - m)!(/ + m')!(/ - m')!] ^ 

- cos (1 + cos (7) 



E 

A;=fc> 



— m — k)\{l — m' — k)\{m + m' + k)\ 



with /i;< := min(/ — m,l — m'), A;> := max(0, — m — m'). The (^1^21™-'!) are the fundamental 
matrix elements, which depend only on the distance R between the two atoms considered; they 
give the matrix elements {hm'l H\l2m') when the connecting vector is aligned with the z-axis, 
R = Re^. 

In [2] the approach of [1] has been extended to 1** and 2"'^ order derivatives of the matrix element 
with respect to the Cartesian coordinates X, Y, Z of the connecting vector R. The approach 
detailed in [2] can be applied to obtain all first and second partial derivatives for all a, (5, R, 
with the exception of /dXdY at the poles f3 = 0, vr. This is due to a coordinate singularity 
and for the specific derivative mentioned no way around this problem was found within the 
approach of [2]; a recourse to a Cartesian description was necessary. The calculations, however, 
were not detailed in that publication. This will be done in what follows. 

The general expression for a Slater-Koster matrix element is a linear combination of i?-dependent 
functions, the (/i^2|"T'|), where the coefficients are combinations of integer powers of the direction 
cosines X/R, Y/R, Z/ R of the connecting vector. Therefore, the general form of a matrix element 
can be expressed as 

M := (/imil/ams) = Yl X^Y' fijk{R) (8) 

ijk>0 

where fij^ is a combination of fundamental matrix elements {hhlml) with constant coefficients, 
divided by R^~^^~^^. At the poles Z = iti?, X = Y = only certain terms from the sum (jH]) will 
remain in OxOy-M- These are precisely those terms where X^Y^ Z^ = XYZ^, as we now show: 



dxM= 

ijk>0 



iX^-^Y^Z^h,k{R) + ^'^'^'^'|^/*^'^(^) 
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Note that this is also vahd if i = 0. 
dydxM = 



ijk>0 



This is also valid if j = and/or i = 0. As i, j > it follows that for X = Y = 0, Z = zizR the 
2nd, 3rd and 4th term vanish. The first term is non-zero in this case only if i = j = 1. This is 
the proposition. 

2 Intermediate calculations 

It is convenient to rewrite the general expression as 

m'=l 

Bm.Bm,TX^,T^,^] {hl2\m'\) (9) 

with S^^^, = AmSl^^, and T^^, = B^T^^,. We have for |mi| > and jmsl > 0: 

-4^1-4^2 = [T(mi)cos(|mi|a) + r(-mi) sin(|mi|a)] x 

[r(m2) cos(|m2|a) + t(— 7712) sin(|m2|a)] = 
^_2'jmi+m2 ]^T-{mi)T{m2) cos(|mi|a) cos(|m2|a) + r(— mi)r(— 7712) sin(|mi|a) sin(|m2|a)-|- 

r(— mi)r(m2) sin(|mi|a) cos(|m2|a) + r(mi)r(— 7712) cos(|mi|a) sin(|m2|a)] (10) 

and 

Bm^Bm2 = [r(-772i)cos(|r7ii|a) -r(mi)sin(|77T,i|a)] x 

[t(— 7712) cos(|m2|a) — r(m2) sin(|m2|a)] = 
( — 1)™'^"'"™^ [r(— mi)r(— 7712) cos(|mi|Q;) cos(|r7i2|a) + r(mi)r(m2) sin(|?TT,i|a) sin(|m2|a) — 

r(— mi)r(m2) cos(|mi|a) sin(|77i2|a) — r(mi)r(— 7722) sin(|r?7i|a) cos(|r?72|a)] (11) 

For 777 7^ we obtain 

AqAjti = — 1=— [r(77i) cos(|777|a) + T{—m) sin(|7n,|Q;)] (12) 
v2 

and BoBm = 0. 

A prefactor XYZ'' in the terms of Q requires an a-dependence that takes the form sin(2a) = 
2 sin a cos a. If we convert the products of two functions cos, sin in H10() . (|ll() into sums of two 
such functions according to the usual trigonometric relations, we obtain sin functions only for 
products of the form sin cos, according two 

sin(|77ii|a) cos(|r772|a) = - sin[(|m,2| + |777i|)a] sin[(|7772| — |r77i|)a] 

1 1 (13) 

sin(|r772|a) cos(|777i|a) = - sin[(|m,2| + |777i|)a] + - sin[(|7772| — |r77i|)a] 
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Henceforth we shall take \m2\ > l^n-il; this is no loss of generality, due to the symmetry prop- 
erties of the inner product. From (|13|) and ()12|) we see that we can get the desired sin(2a) for 

Case A |mi| = \m2\ = 1 

Case B.l |m2| — \mi\ = 2, rrii = 

Case B.2 \m2\ — \mi\ = 2, \mi\ > 
The distinction within case B is warranted for by the distinction in the coefficients Am , Bm for 
m = and m > 0. The r functions in (|lfl)) . (|ll|) . (|12j) further restrict the possible choices of mi 
and m2. In case A, we must have mim2 = —1, and we can choose, without loss of generality 
and without violating previous conventions mi = l,m2 = —1. For B.l we must have m2 = —2, 
and for B.2 mim2 < 0. The matrix elements that can give contributions to the mixed partial 
derivative at the poles thus are: {li,l\ H \l2, —1) , {li,0\ H \l2, —2) , and {li,mi\ H \l2,m2) with 
\mi\ > 0, \m2\ = \mi\ + 2, and mim2 < 0. 

In order to identify these contributions we note that XY = i2^(sin cos a sin a. We can find the 
required derivatives if we extract from Q the factor (sin cos q sin q where possible, evaluate 
the 'remainder' at the poles and divide this result by R^. This is equivalent to dividing by XY, 
and for a dependence XYZ'^, the only one giving rise to nonvanishing contributions, is also the 
same as taking the derivative with respect to X and Y at the poles. The /3-dependence in the 
matrix elements stems from products of Wigner d-functions, as can be seen from Q,®,®- 
The following relations will be needed: 



The d-functions d^^/ occuring in the matrix element all have m > 0, and with p4j) we can 
reformulate all expressions such that m' > and m < m' . We therefore choose to rewrite 
for m, n > as 



dL(m+n)i^) = 2-\-iy-^-- [{I + m)\{l -m)\{l + m + n)\{l -m- n)!] ^ 

< -If {I - cosi?)'-'=-"-t(l + cosi?)'=+'"+5 



X 



E 



^ ^ k\{l — m — k)\{l — m — n — k)\{2m + n + k)\ 



2-'(_l)'-™-" [(/ + m)\{l - m)\{l + m + n)\{l - m - ny]^ 

X (1 - cos^)tfl + cos^)t ' '(1+^°^^)'^"^ . 

x(l cos^yjni + cos^j^ ^,_^i_^_^y_^i_^_^_j.y_^2m + n + k)\' 

(sint?)"(l + cosi?r2-'(-l)'~'"-" [(/ + m)\{l -m)\{l + m + n)\{l -m- n)'l ^ 



' ™ " (-l)'=(l-cosi?)'-™-"-^'(l + cosi?)'= 
^— ^ k\{l — m — k)\{l — m — n — k)\{2m + n-\- k)\ 



(15) 



Only the behaviour of these functions at = 0, vr is important. At "i? = only the term with 
k = I — m — n contributes, and gives (where we don't evaluate sin-i? and cos"!? outside the sum) 



At ?? = TT the only contribution comes from A: = and reads 
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We now look separately at the cases A, B.l, B.2, extracting (sin cos a sin a. 



2.1 Case A : mi = 1, m2 = — 1 

The relevant terms are 

2AiA-id!'iQd!iQ{lil20) = 2 cos a sin ad!'iQd!'iQ{lil20) = 2cosasinadQ\(iQ\(/i^20) (18) 

and 



sinacosa [^t'S^rn' - '^L''^-irn'] ihhlm'l) = 

m'=l 

2sinacosa Yl ("1)"' fc' + (^i^2|m'|) 



m'=l 

l< 



(20) 



2sinacosa [i-^t^'^'^'d'L'W'L'i^ ' P) + i'^)''^^^ ihh\m'\) 

m'=l 

(19) 

In the last equation the following relations have been used 

'^lm"^-lm' ~ "lm'"lm' "l-m'"l-m' ^1 "lm'"l-m' V ^) "l-m'"lm' 

^Im'-'-lm' ~ "lm'"lm' "l-m'"l-m' ^) "lm'"l-m' I ^) "l-m'"lm' 

and ((Tl)) has been applied. 

Using (dH) we get at /3 = for (|T8|l 

cosQsina(sin/3)2^y(I7Tl)Ml7+l)^(Zi^20) (21) 
and (— l)'i+'2x this result at /3 = vr. 

For (|19j) we see from (fT^ that in the case n > the product d!'im'{P)d^im'i'^ ~ l^) contains at 
least a factor (sin/?)^. As we only extract (sin/3)^, the remainder vanishes at /? = and /3 = vr. 
So only the case m' = 1 can produce a contributing term; it is 

sin a cos Q (sin ^) 2 ( - ^[hih + l)+h{h + ^)]{hhl) (22) 

at /9 = and (— l)'i+'2 x this expression at /? = vr. 

2.2 Case B.l : mi = 0, ma = -2 

The relevant terms are 

2^0^-2^0040(^1^20) = 2\/2sinacosado|)d[f2(^i^20) (23) 
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and if /i > also 

l< l< 
AoA^, 4'm'^%m'ihl2\m'\) = ^/2sinacosa Yl ^L'^-2m'ihh\m'\) = 

m'=l m'=l 

l< 

V2sinacosa + <-m'(/3)4-m'(/?)+ 

m'=l 

l< 

V2smacosa ^ + (-iyi+'2d[,;„,(vr - /3)dt,(7r - 

m'=l 

{-ir'^''4UP)dL'i^ - /?) + - (/i/2|m'|) (24) 

Equation gives at /? = 

1 



sinacosa(sin/?)2— ^^(/2 + 2) (/a + iMh - l)(/i^20) (25) 



h+i 



2 



and at /3 = vr we obtain this expression x(— 1)^ 
Turning to (|^H) we first note from (|T^ that for m' > 3 the products d^^^/d^^, contain at least a 
factor (sin /3)^ and thus after extraction of (sin /3)^ still vanish at /3 = 0, vr. This leaves us with 
m' = 1,2. We find: 

4(/?)4(/3) = -d'o\{f3)d[%{P) = -{smPf^^{h + l)h{l2 + 2){h-l) (26) 

at /3 = 0, whereas this product vanishes at /3 = vr after extraction of (sin/3)^ due to a factor 
(1 + cos /?). (iQ\(7r — P)d!'2i{iT — P) vanishes at /3 = due to a factor (1 — cos (3), and at /3 = vr 
results in (|26|) . 

(iQ\(/3)(i2i('/r — (3) vanishes at /? = and at /? = vr gives 

4M4\{7r -P) = -<4(7r - /?) = (-l)'Hsin/3)2iv'(li + 1)^1(^2 + 2)(l2 - 1) (27) 

dQ\(7r — f3)d!'2i{(3) yields p7jl at /3 = and vanishes at /? = vr. 
So from m' = 1 we get a contribution which at /? = is 



sinacosa(sin/3)2 — v/(/i + l)l^{l2 + 2)(/2 - l)(/i^2l) (28) 
V2 

and at /? = vr just (— l)'i+'2 x this result. 
At /? = 

4(/3)42(/3) = i(sin/3)V(Zi + 2)(^i + l)/i(/i-l) (29) 

whereas this product does not contribute at /3 = vr; there ~ f^)d%{'^ — P) gives the above 

expression and in turn vanishes at /? = 0. At /? = dQ2(/?)d22('''" ~ /^) vanishes and at /3 = vr 
yields 

(_l)'il(sin/3)V(/i + 2)(/i + l)/i(/i - 1) (30) 

o 

Vice versa d!'^2i'^ ~ P)d''i2{P) gives the above expression at /? = and vanishes at /? = vr. Therefore 
the m! = 2 contribution is 

sinacosa(sin/?)2^A/(Zi + 2){h + l)li{h - l){hh2) (31) 
2\/ 2 



and at /3 = vr just (— l)'i+'2x this result. 
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2.3 Case B.2 : \mi\ > 0, jmal = \mi\ + 2 

There will be no contribution from the (/1/2O) term in Q, because 

do]m/oU\ (sm/3)l'"^l+l'"^l = (sin/3)2|™il+2 (32) 

and, because \mi\ > 0, the remainder after extraction of (sin/3)^ vanishes at 13 = 0, vr. In the 
sum over (Zi/2|"T''|) for m' > only those parts of the prefactors will contribute to the derivative 
we seek that depend on a like cos a sin a. If we drop all other terms, this 'reduced' sum 
reads in the case mi > 0, 1712 < 

i< 

^ = (-1)'"^+'"^ sin a cos a Yl fen^^^^n^' + (33) 

m'=l 

and in the case mi < 0, m2 > it is just the negative of (jHSI)- We will be using mi > 0, m2 < 
until stated otherwise. With ()14() we can write 

{-iy^ + \m2\+m' Ji ,(vr-/?) + 

\ ' \mi\m' ' \m2\m' ^ ' 

\Zi+|mi|+m' r«i / 



-mm2\rrA(^) ^ 



'^mim'i^)'^m2m'if^) ~ ^|mi |m' (/^)^|m2 |m' (^) + 

(-I)h+l2 + \m^\ + \m2\di ,(7r-P)d\^ , ,(7r-/3)- 
\ ' \mi\m' ^ ^1 \m2\m' ^- 1^ ' 

k(/3) (35) 

Therefore 

c'l c'2 _|_ n-^i n-h _ 
mim' m2m' ' mim' m2m' 

2d\^ , ,(/3)df2 I ,(/?) + 2(-l)'i+'2^fi I Att-P)^^ I ,(^-/3) (36) 

|mi|m'^'^' |m2|m'\'^' ^ ' \mi\m' ^ ' \m2\m' ^ ' ' \ 

Note that (-l)™i+'"2 = (_i)l"^il+l™2| = (_i)l"i2|-|mi| ^ ^ ^.j^g present case. The second 
term in involves the same functions as the first, only evaluated at tt — /3 instead of /3. In 
order to evaluate^ the first term at /3 = 0, we have to choose •& = (3 in ((TH|) : for /3 = tt, = /3 
has to be referred to (|17|) . and we can see that this vanishes (m > 1). For the second term at 
(3 = 0, the required choice is = vr — /? in ((T7|) . This expression vanishes, as [1 + cos(7r — /J)]™" = 
[1 — cos(/3)]'" ^ as /? — > 0. At /3 = TT the second term is evaluated by plugging ^ = -k — (3 into 
((TB|) . This evidently gives the same result as putting /3 = in ((T^ . Looking at H36() as a whole 
it follows that at /? = vr it results in (— l)'i^'2x its value at /3 = 0, the latter being determined 
entirely by the first term. 

If m' > |m2| = |mi|+2 the product d!'\mi\m''^\m2\m' contains at least a factor (sin/3)^. It therefore 
does not contribute. If m' < |mi| = Imal - 2 then , ,d\'' , , = (-l)l™il+l'"2ld'i J\ 

I 1 I ^1 \mi\m' \m2\m' ^ ' m'\m\\ ■m'\m2\ 

^Always explicitly keeping (sin/3)'^. 
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and this again contains at least (sin/3)^. So we have to consider expHcitly the cases m! = |m2|, 
m' = \mi\ + 1, and m' = \mi\, and only at /? = 0. 
m' = \m2\ : 

dfi „ II |(/3) = 



-(sin/3)VGi + kil + 2)^1 + \mi\ + - \mi\){h - |mi| - 1) (37) 
m' = \mi \ + 1 = \m2\ — 1 : 

'^|mi|(|mi|+l)(/^)^|m2|(|m2|-l)(^) = ""^Imi |(|mi|+l) |+I)|m2| = 

m' = \mi\ : 



- (sin/3)2- VGi + |mi| + - |mi|)(/2 + + 2){l2 - - 1) (38) 



dfi II |(/3)df2 II ,(P) = d\' II |(/3)df2 „ m 

\'mi\\mi\^'^ ' \m2\\mi\^'^ ' |mi||mi| |mi||m2|^' ' 
1 



^(sin/3)V(^2 + \mi\ + 2){l2 + \mi\ + I){l2 - \mi\){l2 - |mi| - 1) (39) 

3 Evaluation of (^i^il H \l2Tn2) at the poles 

With the results of the previous section we now compose expressions for the matrix elements at 
the poles. As explained earlier, dropping cosQ;sinQ;(sin/3)^ from these expressions and dividing 
the remainder by gives the derivative we want. Thus we find at /3 = 0: 
Case A: 



dxdY{ll,+l\H\l2,-l) = -^|^v'^l(/l + l)/2(/2 + l) (/1/2O) 



~[/l(Zi + l)+Z2(^2 + l)](/l/2l)}> (40) 



Case B.l: 



dxdy {ll,0\ H \l2, -2) = ^-^V(^2 + 2)(/2 + 1)^2(^2-1) (^1^20) 

- (1 - 5i,o)^^Vhih + m2+2){l2-l) ihhl) 

1 



+ (1 - 6i,o)il - Si,i)^^^^ih + 2)ih + iMh - 1) (/i^22) (41) 



Case B.2 (mi > 0, m2 < 0, otherwise multiply result below by —1) : 



dxdy {h,mi\ H \l2,m2) 
I 1 



^ ^^/il2 + + 2)(Z2 + \mi\ + I){l2 - |mi|)(/2 - - l)(/iZ2|mi|) 



- (1 - Si^\m,\)lVih + + l)(/i - |mi|)(/2 + |mi| + 2)(/2 - |mi| - l)(/iZ2(|mi| + 1)) 
+(l-%|mi|)(l-^/i(|mi|+i))7 VGi + \mi\ + 2){h + \mi\ + l)ih - \mi\){h - - I)ihl2\m2\) 



(42) 
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At /3 = vr we get (— l)'i+'2 x these values. 

Expressions (|4U() and (|41l) were stated in [2]. In this work their derivation has been demonstrated. 
Result (j42j) was not mentioned in [2]. 
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